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Abstract 

We prove that "first singularities" in the non-trapped region of the 
maximal development of spherically symmetric asymptotically flat data 
for the Einstein- Vlasov system must necessarily emanate from the center. 
The notion of "first" depends only on the causal structure and can be 
described in the language of terminal indecomposable pasts (TIPs). This 
result suggests a local approach to proving weak cosmic censorship for 
this system. It can also be used to give the first proof of the formation 
of black holes by the collapse of collisionless matter from regular initial 
configurations. 

1 Introduction 

A fundamental problem in mathematical relativity is to resolve the so-called 
weak cosmic censorship conjecture, the statement that for "reasonable" Einstein- 
matter systems, generic asymptotically flat data do not lead to singularities 
visible from infinity. 

The notion of "reasonable" above is of course not a precise one, and depends 
very much on the context one has in mind. A natural matter source for models 
is provided by kinetic theory. The simplest example is then a self-gravitating 
collisionless gas. The study of the equations describing such a gas, the Einstein- 
Vlasov system, was initiated by Choquet-Bruhat in where the existence of 
a unique maximal development was proven for the Cauchy problem. 

The problem of weak cosmic censorship concerns the global behaviour of the 
maximal development for asymptotically flat initial data. Given the current 
state of the art in nonlinear evolution equations, symmetry must be imposed 
on initial data for there to be any hope of making progress. The global study 
of the initial value problem for the Einstein- Vlasov equations for spherically 
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symmetric asymptotically flat initial data was begun in [7j, where, in particular, 
it was proven that for sufficiently small initial data, the maximal development 
was future causally geodesically complete. The analysis took place in so-called 
Schwarzschild coordinates. In [5], an extension principle was proven, again 
in these coordinates, saying in particular that if the solution stopped existing 
after finite coordinate time t, there was necessarily a singularity at the center. 
These results were meant to provide a first step for a global existence theorem 
in Schwarzschild coordinates. If this coordinate system could then be shown to 
cover the domain of outer communications, and if null infinity could moreover be 
shown to be complete, this would then imply a proof of weak cosmic censorship 
for this system. 

There is another approach to the problem of weak cosmic censorship, due 
to Christodoulou [3j , for the problem of a self-gravitating spherically symmetric 
scalar field. Christodoulou showed that initial data leading to a naked singu- 
larity was codimension 1 in the space of all initial data. This was shown by 
embedding such exceptional data in a one-dimensional subset of the space of 
initial data, such that all other initial data in this subset evolved to a spacetimc 
with the following property, which can be expressed in the language of causal 
sets Given a terminal indecomposable past (TIP) with compact intersec- 
tion with the Cauchy surface, then the domain of dependence of any open set 
containing this intersection contains a trapped surface. The statement that this 
latter property is true for generic initial data can be termed the trapped sur- 
face conjecture. From this property, the completeness of null infinity was then 
inferred, proving weak cosmic censorship. 

It turns out that the relation between the existence of trapped surfaces and 
the completeness of null infinity is quite general. Specifically, in it was 
proven that a weaker version of the trapped surface conjecture is sufficient to 
prove weak cosmic censorship for a wide variety of matter in spherical symmetry. 
In particular, the completeness of null infinity follows from the existence of 
a single trapped or marginally trapped surface in the maximal development. 
The only really restrictive hypothesis on the matter is that "first" singularities 
necessarily emanate from the center. Here, the notion of "first" is tied to the 
causal structure and can be formulated in terms of TIPs. 

The goal of this paper is to prove that the above mentioned hypothesis of 
is indeed satisfied by the Einstein- Vlasov system. As noted before, extension 
principles similar in spirit to this one have been proven before (cf. [EIE3])- These 
earlier results, however, concern the portion of the development of the Einstein- 
Vlasov system covered by particular coordinate systems. Thus, these previous 
results, as far as they concern the maximal development itself, are weaker than 
the results presented here, and in particular, are not sufficient to deduce the 
assumptions of [T2]. 1 

Finally, we make the following remark: In view of UJ, there do exist spheri- 
cally symmetric asymptotically flat initial data for the Einstein- Vlasov system 

^f course, the results of IHIllUI also say something about the behaviour of the coordinate 
system to which they apply, something not addressed here. 
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possessing a trapped surface. Thus, the results of this paper provide in particular 
the first proof of the existence of solutions for collisionless matter representing 
the formation of a black hole. 

2 Initial data 

Initial data in this paper are always given as follows: 

1. We have a C°° Riemannian manifold (E,g), together with an additional 
symmetric 2-tensor K ab , such that there do not exist closed antitrapped 
surfaces in the data, and a compactly supported function fo defined on 
the tangent bundle of E, such that these satisfy 

R - K ab K ab + {tvKf = 16tt / f ( P a )p a Pa/ (1 + p a Pa) 1/2 Vgdp 1 dp 2 dp 3 

V a K a b -V b (tvK) = 8ir J Mp^PaVSdp 1 dp 2 dp 3 

Here the metric g is used to move indices and to define the trace and 
covariant derivative. R is the scalar curvature of g and y/g the square 
root of its determinant. 

2. A smooth S*0(3) action on E such that g, K ab , fo are preserved, and such 
that ~E/SO(3) inherits naturally the structure of a 1-dimensional manifold. 

Here and throughout this paper physical units are chosen so that the grav- 
itational constant has the numerical value unity. We recall the definition of a 
closed antitrapped surface. Let S be a surface in E which is closed, i.e. compact 
without boundary. Suppose that there is a preferred choice n a of an outward 
normal to this surface and let a ab be the second fundamental form of S in 
E corresponding to the outward normal. Then S is said to be antitrapped if 
trcr < -tiK + K ab n a n b . 

3 The maximal development 

The theorem of Choquct-Bruhat pQ, applied to the data considered here, to- 
gether with a standard argument on preservation of symmetry, yields 

Proposition 1. There exists a unique C°° collection (M,g,f) such that 

1. g and f satisfy the Einstein- Vlasov equations 

2. (M., g) is globally hyperbolic, 

3. (A4,g,f) induces the initial data (E, g, K, fo) and E is a Cauchy surface 

4- Any other collection (M,g,f) with these properties 1-3 can be embedded 
in the given one. 
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Moreover, SO(3) acts smoothly by isometry on M. and preserves f, and Q = 
M./ SO(3) inherits the structure of a time-oriented 2-dimensional Lorentzian 
manifold, with timelike boundary T, the center. 

Let 7r : M. — ► Q denote the natural projection. On Q we can define the 
so-called area-radius function 

r(p) = \J Area (it _ 1 (p) ) /47r 

We have r(p) = iff p £ T. We can always choose global future directed null 
coordinates on Q, i.e. such that the metric takes the from —Q 2 dudv. The metric 
of M. then takes the form: 

-n 2 dudv + r 2 7 (1) 

where 7 = jABdx A dx B is the standard metric on S 2 and x A , A = 2,3, are local 
coordinates on S 2 . Let u and v be chosen so that ^ points "inwards" and 
Jjj "outwards" . Such definitions are meaningful in view of the assumption of 
asymptotic flatness. We define 

v = d u r 
X = d v r 

The assumption of no antitrapped surfaces initially means by definition that 

v < (2) 

holds on the initial hypersurface. It follows that it holds throughout Q as a 
consequence of the Einstein equations and the dominant energy condition [2]. 

We shall call the region where A > 0, the regular region, and denote it 1Z. 
We call the region where A = the marginally trapped region, and denote it 
by A, and finally, we shall can the region where A < the trapped region, and 
denote it by T. 

4 The extension theorem 

The extension principle proven in this paper will apply to a region T> C Q with 
Penrose diagram: 




(i.e. a subset V = [111,1*2] x [wi,^] \ (1*2, V2)) such that 

V c K U A. 
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Let Ci n and C ou t be the parts of the boundary of V defined by v = V\ and u = u\ 
respectively. One can think of T> as the "top" of a non-trapped non-central 
indecomposable past (IP) corresponding to a candidate "first" singularity. In 
this language, the result of this paper is that such an IP cannot be a TIP. i.e. 

Theorem 1. IfV C Q, then V C J~(q) for a q G Q. 

The theorem thus says that there is no singularity of this form after all! 

As one might expect, the proof of Theorem^proceeds by obtaining a priori 
estimates in T> and then applying an appropriate local existence result. The a 
priori estimates make use of a certain energy flux along null hypcrsurfaces. This 
fact, together with the fact that regular null coordinates can always be chosen, 
makes it natural to stick to these. We give the form of the equations in local 
null coordinates in the next two sections. Then, in Section [7] we formulate a 
local existence theorem (Proposition |2J) for a double characteristic initial value 
problem. The "time" of existence, in the sense of null coordinates, will depend 
only on the C 2 norm of the metric and the C 1 norm (and the support) of /. We 
obtain energy estimates in Section|Sl and use these, together with the structure 
of the Vlasov equation, to derive in Sections I51 I1UI a priori estimates for the norm 
of Proposition^ The proof of Theoremnwill follow immediately in ScctionlTTl 
Finally, in Section^] we state two applications of our results, discussed already 
in the Introduction. 

The above theorem depends on having a well-behaved matter model and the 
analogous result must be expected to fail for dust. This is illustrated by the 
Penrose diagram Fig. 1 in ^2] ■ 



5 The Einstein equations in null coordinates 

The reader should consult for general facts about the initial value problem 
in spherical symmetry When specialized to this case, the Einstein equations 
are: 

ft 2 1 

d u d v r = — - \v + AirrT UVl (3) 

4r r 

VL 2 1 7rf2 2 
d u d v log ft - -4irT uv + ^ + ^ ~ —1 AB Tab, (4) 

d v (Q-' 2 d v r) = -47rrT OT ,ft- 2 , (5) 



d u {Q- 2 d u r) = -4irrT uu n- 2 . (6) 

The former two equations can be viewed as wave equations for r and ft, while the 
latter two equations can be viewed as constraint equations on null hypcrsurfaces. 
A specific choice of matter model, such as a collisionless gas, leads to expressions 
for the components of the energy-momentum tensor. 
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6 The Vlasov equation 



To describe the Vlasov equation in local coordinates, we need a coordinate 
system on TM. Let p u , p v , and p A denote the functions on TM, defined by 
writing an arbitrary X G TM as 

v „ ^ „ d A d 
X = p u —+p v —+p A 



du dv dx A 

Together with the pull-back of the coordinates on spacetime these functions 
define a local coordinate system on TM. 
Let P C TM be defined by 

P = { 9 (X,X) = -1}, 

where X ranges over future-pointing vectors. We call P the mass shell. It 
follows that 

-n 2 P y v + r 2 1ABP A P B = -i (7) 

We use p u , p A and the pull-back of the coordinates on spacetime to define 
coordinates on P and p v is regarded as a function of these coordinates defined 
by the relation Q. The Vlasov equation is an equation for a non-negative 
function 

which, in the case that / is spherically symmetric, is given by 

p u ^+p v % = {d u (\ og n 2 ){p-f + 2n-\\ 1AB p A p B )^L 



-Xp v )p a ^a- 



In deriving this we have used the expressions for the Christoffel symbols given in 
Appendix A and the fact that a spherically symmetric function / on the mass 
shell is a function of the variables u, v, p u , and r )ABP A p B - This implies the 
identity 

p dx A ~ LbcP p 8p A 

which has been used to simplify the Vlasov equation. Note that both the ex- 
pressions "fABP A P B and p A -g^x have a meaning independent of the particular 
choice of coordinates x A on S 2 . 

Finally, to close the system, we must define the energy-momentum tensor. 
We first note that for any point q £ M, it follows that P q , as a spacelike hy- 
persurface in T q M, inherits a volume form from the Lorentzian metric. In local 
coordinates this volume form can be written r 2 (p") -1 dp u y/ydp A dp B or alterna- 
tively r 2 (p v )~ 1 dp v ^/jdp A dp B , where y^y is the square root of the determinant 
of jab ■ We then have 

r 2 p aPb f(p u r 1 Vjdp u dp A dp B , (9) 
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where p a = g a bP b ■ It follows immediately that this matter model satisfies the 
energy conditions: 

T uv >0,T vv >0,T uu >0 (10) 

7 A local existence theorem 

To prove our extension theorem, we will certainly need to appeal to some sort 
of local existence theorem. In particular, it is the norm in this theorem that will 
tell us what quantities we must bound a priori in T>. In principle, one could try 
to prove estimates so as to apply the local existence result of For various 
reasons, however, the following local existence theorem for a characteristic initial 
value problem will be more convenient: 

Proposition 2. Let k > 2. Let f2, r be positive C k -functions defined on [0, d] x 
{0} U {0} x [0, d], and let f be a non-negative C k ~ 1 function defined on the part 
of the mass shell over [0,d] x {0}U {0} x [0,d]. Suppose that equations © 
hold on {0} x [0, d] and [0, d] x {0} respectively, where T uu and T vv are defined by 
(01, and suppose in addition that the C k compatibility condition holds at (0,0). 
Define the norm: 

N u = sup {IQUO-^l^fiU^OUrUr-rM^rU^rl, 

[0,d]x{0} 

5,1/1,1^/1,1^/1,1^/^}, 

n v = sup {|n|,|n- 1 |,|a n|,|aSn|,H,|r|- 1 ,|a r|,|^r|, 

{0}x[0,d] 

s.I/USb/I.I^/I.IV/It}. 

N = sup{N u ,N v }, 

were S denotes the supremum of (p u ) 2 + {p v ) 2 + JabP A P B on the support of f 
and \va\j = (i^vaub) 1 ^ 2 ■ Then there exists a 8, depending only on N, and 
C k functions (unique among C 2 functions) r, f2 and a C k ~ 1 function (unique 
among C 1 functions) f, satisfying equations (01, Q, (@, ©, © in [0,(5*] x 
[0,(5*] , where S* = min{(i, S}, such that the restriction of these functions to 
[0, d] x {0} U {0} x [0, d] is as prescribed. 

Proof. See Appendix iBl 

The compatibility conditions referred to in the statement of the proposition 
are as follows. The data includes the values of the function / on the part of the 
mass shell over [0, d] x {0}. All derivatives of / tangential to this manifold can 
be calculated by direct differentiation. By using the field equations transverse 
derivatives (and thus all derivatives) of / can be computed up to order k— 1. In 
a similar way, all derivatives up to order k — 1 can be computed on {0} x [0, d]. 
The condition that derivatives determined in these two different ways agree at 
(0,0) is what is referred to above as the C k compatibility condition. 

Let us add the remark that, defining g on M. by Q, the above gives rise to 
a solution of the Einstein- Vlasov equations upstairs, with the obvious relation 
to characteristic data, interpreted upstairs. 
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8 Energy estimates 



A fundamental fact about the analysis of spherically symmetric Einstein matter 
systems in the non-trapped region is the existence of energy estimates. 

To describe these, let us first settle for a particular null-coordinate descrip- 
tion of the set T>. We normalize our u-coordinate such that v = —1 along Cj n . 
For the v coordinate, we first define the quantity 



1 

4 



= --n 2 -- 1 



and then define v such that k = 1 along C ou t- T> is thus given by [0, U] x [0, V] \ 

mv)}. 

The concept of energy in spherical symmetry is given by the so-called Hawk- 
ing mass, given by: 

m= 7 -(l- 8 a rd a r) = T -{1 - 2g uv d u rd v r) = ^(1 + 4fT 2 A^). 

We will also introduce the so-called mass-aspect function 

_ 2m 
^ r 

Note that 

k(1-h)=X. (11) 
From l|3Jl 10, we compute the identities: 

d u m = %irr 2 Vl~ 2 (T uv v - T uu \) 

= -27r K -V 2 T M , + 2Tr^—^r 2 T uu (12) 



d v m = 87rr 2 fi 2 {T UV \ — T wu v) 

= -2ir^—^r 2 T uv + 2iTK- 1 r 2 T vv . (13) 



The first point to note is that the signs of (|12fl and (|13f) . together with the 
signs of A and give a priori bounds for both r and m. Indeed, set 

m = m(U, 0) > 0, 

r = r{U, 0) > 0, 
M = m(0,V), 
R = r(0,V). 

By © and the fact that V C TZ U A, we have that 

ro < r < R (14) 
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throughout V. On the other hand, l(T2|l . (|T3|l and ((Tn|) give d u m < 0, > 0, 
and thus 

mo < m < M. (15) 

Now we make a trivial observation. In view of the fact that we have the a 
priori bounds JTSJ, if we reexamine the equations (|12fl . I|13|) . keeping in mind 
that both terms on the right hand side have the same sign, we obtain the bounds: 

2?r ( 1 Z V) r 2 T uv (u, v)dv < M — m , (16) 

r-2 



2TTK~ 1 r 2 T vv (u, v)dv < M - m , (17) 

ui 

112 



/ 2iTK~ 1 r 2 T uv {u,v)du < M - m , (18) 

J Ui 

(u,v)du < M - m . (19) 

These will be our energy estimates. 

As we shall see, our use of the above estimates will not quite be symmetric 
for u and v. The reason is this: The "constraint" equation © can be seen to 
be equivalent to the following equation for k: 

3 u k = Anrv- 1 T UU K. (20) 

From J2J, (|20fl and (jl()|) . we see immediately 

0< k < 1 (21) 

throughout T>, i.e. k" 1 > 1. This means that a priori we control jT vv dv, but 
not J T uu du. 

Finally, note that we can rewrite equation © as 

d v v = 2r- 2 nvm + AirrT uv , (22) 

or alternatively 

d u X = 2r~' 2 Kvm + AirrT uv . (23) 
Thus, integrating (JUJ), in view of (J2U, CSJ), (UU), and (JTUJ), we have that 

u > - e 2r o 2MV = -N. (24) 
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9 C estimates for the metric 

So far, we have not used the Vlasov equation, only the energy condition (|10|) . 
Indeed, all estimates obtained so far are familiar from the results of To go 
further, we must use the Vlasov equation itself and the special structure of the 
energy-momentum tensor. In this section, we shall estimate the support of / 
and show C 1 estimates for the metric. 

Before proceeding, let us give names to bounds on certain quantities on the 
initial segments C; n U C out . Define 

G = max < sup |9„logS! 2 |, sup | <9„ log f2 2 1 > 

[[0,C/]X{0} {0}x[0,V] J 

F = sup /, 

^r 1 ({0} x [0,V]U[0,U] x {0}) 

where tt\ denotes the projection from the mass shell, define £ to be supremum 
of the radius of support of / in the p° and p u directions along n^ 1 ({0} x [0, V] U 
[0, U] x {0}), and define X be the supremum of r )ABP A p B over the support of 
/■ 

Let us note first two easy bounds. Clearly, 

< / < F 

throughout the mass shell over T>. Moreover, by l|14[) and conservation of angular 
momentum applied to geodesies, it follows that 

lAB p A p B {x) < X (25) 

for any x G P in the support of / over T>. In particular, in the expressions 
defining energy-momentum, we can thus always replace an integral over the 
variables p A by the integral over the ball of radius X about the origin. 
We have the following: 

Lemma 1. The inequality 

—Quv9 AB Tab < 2T U1J 

holds throughout T>. 

Proof. The inequality is equivalent to the statement that the trace of the energy- 
momentum tensor is non-positive. This holds for collisionless matter indepen- 
dently of symmetry assumptions. It is proved straightforwardly by taking a 
trace in the formula defining the energy-momentum tensor in general coordi- 
nates with the spacetime metric. □ 
We can rewrite as 

d u (d v \ogfl 2 ) = -87rT uv -4Kmr- 3 v + 8nKvr- 2 <y AB T A B- (26) 
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Integrating H26|) . applying the above lemma, the energy estimate (|18|l . and the 
bounds dJ, (ESJ, (|2J, we estimate <9„logf2 2 : 



|<9„logft 2 | < G 



SttT uv cIu — / Sirnvr 2 "/ ab TabcIu 



Anmr vdu 



< G+ / Anr 2 {2nr 2 n 1 T uv )du- / Anmr^vdu 



< G + 4r 2 J 27rr 2 K^T^du - J 4 

< G + 4r 2 y 27rr 2 K- 1 T tll ,dw - / 4 



< G + Ar 2 (M -m Q ) + 2(r^ 2 -R-' 2 )M 

= G'. (27) 

Integrating now l|42|l. using l(2*7|) . we obtain 

|logft 2 ( w , v )| < |logft 2 (u,0)|+G'F, 

and thus, since | log0 2 (u, 0)| < C for some C, we have, 

< c< n 2 (u,v) < D. (28) 

Now, we turn to estimate the projection to the p"-axis of the support of 
/. We proceed by considering the geodesic equation. Let -f(s) be a geodesic 
crossing {0} x [0, V] U [0, U] x {0} at s = 0, such that 7'(0) is in the support of 
/. Let p v (s) denote the ^ component of the tangent vector of 7. We have 

W)\s) = -Tl v (p*f-T AB p A p B . (29) 

using the Christoffel symbols in Appendix 1X1 Integrating 129|) . we have now by 



'(s) = p v (0)e-^KAp v )ds _ r v AB p A (S)p B (S)e-^ sr ^^ds 

Jo 

j/ n \ - f *,S r* ito f rt A B - f "/- s ? r" ,du ,~ 
(0)e Jl, (°) — / 1 ^4 B p p e Jl, < s ' as 
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«(0) e -/"$ r - d,J + f ^ 2(- l /)fi- 2 »- 7AB A B e - / "w r «*'( P 1 ')- 1 dt;. 



/«(o) 

Thus, for s' < s, (replacing with s') we have, by (|28|l and l|24|l . the inequality 

P v (s) < p v ( S ')e- j vt>\ T - dv + f" 2Nc- 1 r lABP A p B e- V ™ dv {p")- 1 dv. 

(30) 

Suppose p"(s) > 2S for some < v(s) < V, and let s' be the last previous time 
s > s' > such that p v (s') > 2S, i.e., we have p u (s*) > 2E on [s*, s]. By l|5U)). 
H27|l. the angular momentum bound (|25fl . and (|42|l . we have 

p*(«) < 2£e yG ' + 2RX 2 Nc- 1 e VG 'v(2T l )- 1 , 
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P v (s) < C. (31) 
We can now easily estimate T uu pointwise: 

T U u = T [ r 2 {p u ) 2 f%^dp A dp B 

JO J\1ABP A P B \<X P 

= (.9«,) 2 r I r 2 (p^f^V7d P A dp B 

JO J\1ABP A P B \<X P 

2. .2 f f „2/„ v\2 fdp" i— , A, B 



4«/V / / fWYf^-yftdp**? 
JO J\ 1AB p A p B \<X P 

\v 2 Fk 2 [ C [ r 2 (p v )dp v ^dp A dp B 

JO J\~iARV A P B \<X 



'0 J\ 1AB p A P B \<X 

< lQnR 2 v 2 FC 2 X 2 = v 2 E 

< N 2 E, 

in view of (|3*T|> . (fS^jl . J5lJ, iffl)) and the angular momentum bound (|2*5|) . (Note 
that T uu v~ 2 < E is a coordinate invariant 2 bound.) Integrating (|2C)|> . we obtain 
now 

K > e ~ J iwr2 i$ Lvdu > e ~^RENU _ 

(Actually, we have in fact already estimated k from below since k^ 1 — A(—v)Vl~ 2 . 
From the inequality 

PuPv < ^{pI+pD, 

we have 

This allows us to estimate 9 U logfi 2 = r^ u : 



iril < G 



8TrT uv dv — I 8irK,vr 2 ~f AB TAB<iv 



< G + AttN'EV + J AirT vv dv - J iKmr~ s vdv 

< C 

We can easily obtain an estimate now for T vv . X can be bounded by integrating 

10 C 2 estimates for the metric 

In this section, we derive C 2 estimates on the metric and C 1 estimates for /. 
The ideas of this section originate in [7]- 

2 i.e. it does not depend on the normalization of u 
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It has already been shown that the following quantities are bounded: r, r -1 , 
m, m , k, k^ 1 , v, v^ 1 , A, f2, fi -1 , all first order derivatives of f2, all components 
of the energy-momentum tensor, and all Christoffcl symbols in 13 7|) 1)42(1 . From 
these estimates and (1221) and (12. ''>t . it follows that d v v and d u X are bounded, 
from l|12|) and (|13[) it follows that d u m and d v m are bounded, and from (|26|) . 
it follows that d u d v ft is bounded. Writing i/ = — -jil 2 ^ 1 and differentiating 
in u, we see from (12011 that d u v is bounded, while writing k = — ifi 2 ^ -1 and 
differentiating in v, we see that d v n is bounded, and thus, from (JTTJ, we see 
that c\A is bounded. These estimates and the formulas Q37 |) -I|42 [l allow us to 
control all first order derivatives of the Christoffcl symbols, except d u Y™ u and 

Q -nv 

Uy-L vv . 

Since the components of the curvature tensor can be expressed in terms of 
those derivatives of the Christoffel symbols which have already been estimated, 
wc obtain bounds for all components of the curvature tensor in our coordinates. 
The above estimates allow us to estimate the first derivatives of the exponential 
map on the tangent bundle. This, in turn allows one to estimate the derivatives 
of / in terms of initial data. 

We can, however, argue more directly as follows. Let us abbreviate the 
Vlasov equation JHJ by X(f) = where X is the Vlasov operator written in 
these coordinates. Note that p v is to be thought of as expressed in terms of p u 
and p A via the mass shell condition (J7J|. Define fi = d u f —p u d u log VL 2 d p ^f. Dif- 
ferentiating the Vlasov equation with respect to v, p u and p A gives the following 
equations: 

x{d v f) = -(d vP v )d v f + (d u d v io g n 2 { P u ) 2 + d v {~2n~ 2 r\) 1ABP A p B )d p vf 

+ 2{d v {vr- l )p u + 9„(Ar- V + \r' l d vP v )p A d pA f (32) 
X(d p uf) = -d u f - (d p »p v )d v f + 2d u logOV^/ 

+ 2{ur- 1 + \r- l d^p v )p A d pA f (33) 

X(p D d p of) = -p D {d p op v )d v f-m- 2 r\ lAB p A p B d pU f (34) 

+ 2r- 1 Xp D d p np v p A d p Af. (35) 

Differentiating the Vlasov equation with respect to u gives the following equation 
for f x : 

X(fi) = -p u d u (logtf)X(d p uf)-d uP v d v f 

+ {- P u p v d u d v log n 2 - d u log n 2 {d u log n 2 { P u ) 2 + 2n- 2 r\ 1AB p A p B ) 

- 2d u (n- 2 r\) lABP A p B )d p ,f 

+ 2(9 tl (^r- 1 K + d u (Ar- V + d uP l 'Xr'^p^A j '. (36) 

The quantity X(d p uf) can be substituted for by one of the previous equations 
and d u f may be eliminated from the equations in favour of j\. The result is 
a linear system of equations for the evolution of {fi,d v f,d p ^f,p d p Af) along 
the characteristics of the Vlasov equation. The coefficients are known to be 
bounded and so we can conclude that d u f, d v f, d pU f and p A d p Af are also 
bounded. (Note that since p u and p v are bounded the derivative with respect to 
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X is uniformly equivalent to a derivative along the characteristic with respect 
to u or v as parameter.) 

From this, we immediately estimate d u T a \> and d v T a t pointwise. We now 
estimate d u T!^ u by differentiating (|26[) in u and integrating in v, and similarly, 
d v T^ v by differentiating in v and integrating in u. Note that |9 p a| 7 can also be 
bounded. This can be seen by passing from polar to Cartesian coordinates and 
noting that the resulting metric components are C 2 . As a consequence / is C 1 . 

11 The Proof of Theorem Q] 

Let N/2 denote the sup of the norm defined in Proposition where the sup is 
taken now in all of T>. By the estimates of the previous section, we have that 
N/2 < oo. Let S be the constant of Proposition[21corrcsponding to N. Consider 
the point (U — S/2,V — 8/2). Translate the coordinates so that this point is 
(0,0). Since Q is by definition open, by continuity, there exists a <5 > <5* > 8/2 
such that 

{0} x [0, 8*] U [0, 8*] x {0} c Q 

and the assumptions of Proposition hold on {0} x [0,8*] U [0,8*] x {0}, with 
N and 8* as already defined. It follows that there exists a unique solution of in 

£ = [0,8*] x [0,8*]. 




Thus the solution coincides in £ n Q by uniqueness. One sees that £ U Q is 
clearly the quotient of a development of initial data. By maximality of M, we 
must have £ U Q C Q. Thus, in particular, in the old coordinates we have 
(Z7, V) G Q, and the theorem holds with q = (U, V). 

12 Applications 

We will say that a spherically symmetric maximal development has a black 
hole, if T + is complete in the sense of 0], 3 and if J~(I + ) has a non-empty 
complement. 

We have shown that the results of apply to our matter model. In par- 
ticular, the fact that the complement of J~(2 + ) is non-empty implies the com- 
pleteness of null infinity. That this set is non-empty can be inferred in turn from 

3 See 1121 for a definition of X+ in this context. 
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the existence of a single trapped or marginally trapped surface. Asymptotically 
flat spherically symmetric solutions of the Einstein- Vlasov system possesing a 
trapped surface were constructed in Thus we have 

Corollary 1. There exist solutions of the Einstein- Vlasov system which develop 
from regular initial data and contain black holes. 

The fundamental open question in gravitational collapse is to show that 
gcncrically, either the solution is future geodesically complete or a black hole 
forms. In view of |12j and the results of this paper we have 

Corollary 2. Suppose that for generic initial data, the maximal development 
either contains a trapped surface or marginally trapped surface, or is future 
causally geodesically complete. Then weak cosmic censorship is true. 

Thus, weak cosmic censorship can be reduced to a slightly weaker version of 
Christodoulou's trapped surfaces conjecture. As remarked in the Introduction, 
this suggests a local approach to its proof (cf. PJ). 
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A The Christoffel symbols 



Note: 



guv 




g uv = -2rr 2 
n 2 = -akv. 



The nonvanishing Christoffel symbols are given by: 



r 



= -g uv r\~{AB 



(37) 



n B = -g m nnAB, 



(38) 




(39) 



1 Bu — vr °Bi 



(40) 



r 



cUogft 2 



(41) 
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C=cUogfi 2 . (42) 

In fact the Christoffel symbols T^ B , which depend on a choice of coordinates 
on the spheres of symmetry need not vanish but the expressions for them are 
not needed in this paper. 

B Proof of Proposition [2] 

The proof of local existence follows from simpler considerations than the proof 
of the estimates of Sections I8HTU1 In particular, one does not need to consider 
energy estimates, for one can recover naive pointwise estimates using the small- 
ness parameter. As in Section^] the idea of [7] again makes its appearance, to 
show C 1 bounds on / directly from C° bounds on the curvature, before bound- 
ing the C 2 norm of the metric. Since all these methods have appeared before, 
we will only sketch the details here. 

Let initial data be fixed. Define the space 

Ac C 2 ([0,S} x [0,8]) x C^fltyi] x [0,(5]), 

for S to be determined later, consisting of all twice continuously diffcrentiablc 
nonnegativc functions r, continuously diffcrentiablc nonncgativc functions fl, 
extending the prescribed values, such that 

AT 1 / 2 < r < 2N, (43) 



N~ 1 /2 < < 2N, 



S up{\d u r\,\d v r\,\d 2 u r\,\d 2 v r\}<2N 



(44) 
(45) 



sup{|5 u fi|, \d v n\} < 2N. (46) 

Consider the subset B C A, consisting of those (r, fi) for which fl is C 2 , and 
for which 

snp{\d 2 u n\, |flgfi|, \d u d v Q\} < 2N. (47) 

Note that the closure of B in A, denoted B, consists of (r, f2) such that d u Q, 
<9„f2, are Lipschitz, with Lipschitz constants given by the above. 

We shall define in the next few paragraphs a continuous map 4> : B — > A 
taking (r, 0) to (f, fi). 

Given r, CI, first, let / be defined to solve the Vlasov equations on the metric 
defined by r and fl, with given initial conditions. Note that since the Christoffel 
symbols of this metric are Lipschitz, it follows that geodesies can be defined, 
and thus / can be defined by the requirement that it is preserved by geodesic 
motion. It follows immediately that 

0<f<N, (48) 
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and, after appropriately restricting to sufficiently small S, it follows easily by 
integration of the geodesic equations that 



S < 2N. (49) 

In the case where (r, O) G B, we have that / is in fact C 1 , since the exponential 
map is differentiable. If 6 is chosen sufficiently small, it is clear from (|43|) (14 711 
that, in this case, we can arrange for 

snp{\d v f\, \d u f\, \d p J\, |cy /| 7 } < 2N. (50) 

Given now /, we can define T uv , T vv , T uw in the standard way. In view 
of (|4*5|) - (|4T)|> . and (|4*5)l . these terms can be estimated. Now, set v = d u r, 

A = d v r. We define f by 

f u C" 1 1 

f(u, v) = r(u, 0) + r(0, -u) - r(0, 0) + / / — r~ 2 Sl 2 \v + AnrVl i T uv dudv 

Jo Jo 4 r 

(51) 

By appropriate differentiation of 151|) . it is clear from our bounds thus far 
that we can define and estimate v = d u f, X = d v f, and d u d v f. We can retrieve 
the bound (|43|l for r by integration of the after restricting to small 6. For 
(r, O) G -B C B, it is clear we can also define and estimate d^r, d%r, by differen- 
tiating H51f) twice in u or twice in v, in view of the fact that all other derivatives, 
including d u T uv , d u v, etc., are clearly defined and bounded, in view of J5UJ), and 
since these derivatives are defined initially. By appropriate choice of 8, we can 
clearly arrange-for (r, f2) £ B-so as to retrieve the bound H45|) . 

Define now Cl > by the relation 

log ft 2 = log ft 2 (u, 0) + log ft 2 (0,v)- log ft 2 (0,0) (52) 

+ / / (-8wT uv + -n 2 ?-' 2 + 2?~ 2 \v -2nn 2 ?- 2 1 AB T AB )dudv. 
Jo Jo 2 

Again, for small enough 5, it is clear that one can arrange for ft to satisfy (|44|l . 

Differentiating (|52|l appropriately, in view of the initial conditions for f2, it 
follows that, for (r, f2) G £?, 17 is C 1 , and for 5 small enough satsfies (|46|l . while 
for (r, 17) G 17 is C 2 , and for 5 small enough, satisfies (|4Tjl . 

Thus, we have shown that after judicious choice of 5, <f> maps B to itself. By 
continuity, it maps B to itself. 

The map 4> can easily be shown to be a contraction in B for the norm of A, 
i.e., we can show that 

d A ((f!, Or), (f 2 ,0 2 )) < ed A ((n,Oi), (r 2 , 2 )), (53) 

for an e < 1 and all (n,Cli) E B. To see this, define first fi, corresponding to 
(ri,Oj). Let Ti denote an arbitrary Christoffel symbol for (rj, Oj). We clearly 
have 

|r x -r 2 | <cdx((ri,fii),(r2,n 2 )). 
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We easily obtain 



\fi-f 2 \<C6sup\r 1 -r 2 \ sup (\dfi\ + 

T 8=1,2 

Clearly we can also bound sup |T™ — J 1 ™! < Csup|/i — /2I. One bounds 
(isi — V2) by expressing d v {v\ — V2) as a linear combination of f2i — Q21 t\ — ^2, 
i^i — ^2, Ai — A2 and (T™ — T^") with bounded coefficients. One immediately 
obtains a similar bound for sup \f\ — f 2 |. The terms sup \d u fi — d u r2\, sup \d v f\ — 
d v ?2\, and sup \d u d v fi — (9 M 9„f 2 |, can be handled in the same way. One then 
obtains a bound of the above form for sup \ d v log Cl± — d v log J7| | , and similarly 
for sup \ d u log f7f — d u log £1| | . Either of these bounds of course implies a bound 
for sup — 

To bound sup \d 2 fi — c?^ r 2 1 , we compute 

+ 47ra u (r0 4 )T u ' 1 ' + 47rrfi 4 d u T m 'dv (54) 
= dlf\ v=0 + J d u (- \r- 2 ^ 2 - r^Xu) + ^d u {rtf)T uv 

- 4Trrn 4 d v T vv + Airril 4 ^ T ■ T)dv 
= flgf |„ =0 - 47rrfi 4 T™(M, v) + 4 7 rrf7 4 T m '(u, 0) 

+ Jo ^ (~i r_2fi2_r_lA ^) +^d u {rn i )T uv 

+ ATrd v (rn 4 )T vv + 47rrfi 4 (^ T ■ T)dv. (55) 

Here we have used the equation V a T ab = 0, which follows from the Vlasov 
equation, and we have integrated by parts. It is now clear that estimates 
for differences follow as before. We argue in an entirely analogous way for 
sup _ d%r 2 \. 

After restricting to sufficiently small 8, all constants in the above bounds 
can be made small. We thus have indeed shown l|53|) . It follows by continuity 
that $ is also a contraction on B C A, and thus, since B is closed, has a fixed 
point in B. 

Given such a fixed point (r, S7), define / as before. To show that (r, f2, /) 
corresponds to a solution of the equations, we have basically only to show that 
/ and d u fl, d v fl, which a priori are Lipschitz, are in fact C . (In particular, 
from this it will follow that the constraint equations (01-© ar e also satisfied.) 
But, in view of the fact that / is initially C 1 , it follows that / is C 1 if the 
exponential map is C . (The C 2 compatibility condition is used at the point.) 
But this latter fact follows from the continuity of the curvature, as shown in 
Exercise 6.2 of Chapter V of jH] 4 - That the curvature is continuous follows 

4 If the reader does want to apply to this fact, then one can argue as follows: in view of 
the computations above, in the space B, we have that curvature is in fact C 1 with estimates; 
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by computation, since r is C 2 , ft is C 1 and d u d v Q is C°, and <9 2 f2 and d^tt 
do not appear in the expressions for curvature. From the C 1 property of /, 
the C 2 property of fl follows immediately. Similarly, higher regularity follows 
immediately if it is assumed. □ 
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